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1.  INTRODUCTION  AND  SUMMARY 


The  Polywell**  is  a  magnetic  version  of  a  Spherically  Convergent  Ion  Focus 
(SCIF)  device  which  was  proposed  by  R.  W.  Bussard1*2  as  a  significant  variation 
of  earlier  studies  on  electrostatic  confinement.3,4  The  idea  of  this  device  is 
to  inject  high  energy  electrons  into  a  quasi-spherical  magnetic  field;  the 
electrons  create  a  potential  well  of  sufficient  depth  to  accelerate  ions  from 
low  energy  at  the  periphery  to  fusion  energies  within  a  focus  at  the  center  of 
the  sphere.  Injection  of  electrons  keeps  the  system  electrically  nonneutral . 
so  that  the  potential  well  which  accelerates  the  ions  is  maintained  at  a  constant 
value  sufficient  to  confine  the  ions  within  the  device,  returning  them  again  and 
again  at  high  velocity  to  the  central  focus.  Essential  to  the  success  of  the 
scheme  is  that  the  ions  maintain  their  nonthermal  velocity  distribution,  with 
primarily  radial  flow,  long  enough  to  produce  fusion  in  the  dense  focus  at  the 
center  of  the  sphere. 

If  the  ion  source  is  at  the  periphery  and  ions  return  to  their  birth  point 
(turning  points),  then  ar  isotropic  source  at  R  will  clearly  produce  an 
anisotropic  (non-Maxwellian)  distribution  at  r  <  R,  since  vr  increases  as 
y2e0(r)/m^ ,  while  from  conservation  of  angular  momentum  v±  increases  as  v  • 
vj.R  R^r’  w^th  2e^max^mi  "  vj.R  R  /rc  determining  rc»  the  closest  approach  of  an 
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ion  to  the  center,  where  v±R  is  the  velocity  perpendicular  to  if  of  an  individual 
ion  at  r  -  R,  and  where  i  »  0  at  the  location  of^he  ion  source.1,2 

If  the  ions  must  converge  to  a  high  density,  small  radius  focus,  to  form 
a  reactor  relevant  configuration,  ion-ion  collisional  effects  can  destroy  the 
configuration,  even  if  they  do  not  cause  ions  to  leave  the  device.  There  are 
two  distinct  effects  that  collisions  have  on  ion  convergence.  First,  as  ions 
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born  at  low  energy  pass  through  the  dense  core  of  the  device,  they  will  be  nearly 
monoenergetic,  at  the  energy  of  the  well.  Collisional  processes  will  tend  to 
thermalize  the  velocity  distribution  in  the  core.  An  ion  which  is  upscattered 
in  energy  will  have  a  turning  point  at  a  radius  larger  than  the  radius  at  which 
it  was  born  and  can  eventually  be  lost  to  the  dense  core  by  escape  from  the 
device  or  by  magnetic  deflection,  since  B  increases  with  radius  in  the  multicusp 
field  characteristic  of  SCIF  devices,  B  -  rm,  with  m  typically  >  3. 

The  second  collisional  effect  involves  isotropization  of  the  ion  velocity 
distribution  in  the  bulk  of  lower  density  plasma  outside  the  core.  In  this 
large  scale  region,  ion-ion  collisions  can  transfer  momentum  from  the  primarily 
radial  flow,  which  produces  good  spherical  convergence,  to  a  more  isotropic 
distribution  of  velocities.  This  deflection  by  collisions  leads  to  an  increase 
of  local  azimuthal  velocity  in  this  bulk  region,  which  clearly  decreases  the 
convergence  of  the  ions  to  the  core. 

In  addition  to  the  core  and  the  bulk  plasma  which  surrounds  it,  there  is 
a  third  region  with  distinctive  collisional  properties,  namely  the  edge  plasma, 
where  the  bulk  of  the  ions  have  their  turning  points.  In  this  region  ion 
velocities  are  low,  and  collisions  much  more  frequent  than  in  the  interior.  In 
this  region  the  ion  distribution  due  to  previous  collisions  in  the  interior  tends 
to  become  more  anisotropic,  while  collisions  in  the  edge  itself  make  the  edge 
distribution  more  isotropic.  An  isotropic  distri  *  ion  at  the  edge  will  produce 
exactly  the  strong  radial  convergence  required  for  the  SCIF.  Thus  due  to  the 
large  ion  orbits,  collisional  processes  in  the  edge  will  tend  to  drive  the 
distribution  in  the  interior  to  stay  anisotropic  and  convergent.  In  this,  edge 
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collisions  compete  with  local  collisions  in  the  Interior  which  attempt  to  produce 
a  local  Maxwellian.  If  edge  collisions  dominate,  even  though  large  angle 
scattering  will  eventually  thermal ize  the  entire  system,  the  more  rapid  small 
angle  scattering  (which  competes  with  the  fusion  reacting  time  scale)  will  force 
the  distribution  to  stay  non-Maxwellian  in  the  Interior.  This  effect  is  simply 
a  matter  of  competitive  rates  of  thermal izat ion  in  the  various  regions. 

In  Section  II,  we  calculate  the  ion  collision  frequencies  in  the  various 
spatial  regions,  for  both  experimental  and  reactor  grade  plasmas.  In  Section 
III,  we  use  a  diffusion  approximation  to  give  heuristic  estimates  of  the  time 
soales  for  "loss"  of  ion  convergence  due  to  ion-ion  collisional  effects  in  the 
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core  and  bulk  of  the  device,  as  discussed  above.  In  Section  IV,  we  give  a 
Fokker-Planck-like  calculation  of  the  Ion  distribution  which  results  from  the 
combined  effect  of  collisions  In  all  three  regions,  and  show  that  the  edge 
effects  dominate,  producing  a  Maxwellian  edge  and  a  non-Maxwellian  convergent 
Interior.  In  Section  V  we  apply  the  results  to  estimating  loss  rates  in  a  SCIF 
reactor. 


II.  ESTIMATES  OF  ION  COLLISION  RATES 


In  this  section  we  estimate  the  various  Ion  collision  rates,  as  discussed 
In  Section  I.  Because  collisions  depend  on  energy  and  density,  we  first  present 
the  radial  density  and  energy  profiles  expected  In  an  Idealized  version  of  the 
SCIF,  and  discuss  the  colllslonal  rates  In  the  different  regions,  for  both 
experimental  and  reactor-grade  parameters,  using  the  parameter  set  In  Table  I. 


A.  ION  ENERGY  AND  DENSITY  PROFILES 


As  a  first  approximation  for  the  motion  of  a  test  Ion  In  a  SCIF  device, 
we  assume  the  Ion  undergoes  conservative  motion  In  a  spherically  symmetric 
electrostatic  potential  well.  The  ion  Is  born  with  an  energy  EQ  at  a  location 
R.  The  total  energy  of  the  particle  Is  conserved: 

E^r)  -  \  mi  (v*  +  v*)  -  etf(r)  .  (1) 

Here  E-  is  the  total  energy,  m^  is  the  ion  mass,  vr  and  v±  are  the  radial  and 
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azimuthal  ion  velocities,  v*  ■  v£  +  v^  In  a  spherical  coordinate  system,  and 
d(r)  Is  the  electrostatic  potential  energy. 

We  assume  the  electrostatic  potential  has  the  following  spatial  form: 

-  *«»x  [>  -  <S>P]  '  <2> 
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where  4max  Is  the  value  of  the  potential  at  r  •  0,  the  center  of  the  device,  and 
4  •  0  at  r  -  R,  where  R  Is  the  location  of  the  Ion  source.  The  actual  form  of 
4(r)  Is  not  relevant  to  the  qualitative  results  we  obtain. 

Ue  assume  that  the  Ion  Is  born  at  a  location  where  magnetic  effects  are 
weak.1'2  If  this  Is  not  the  case,  the  effect  of  collisions  will  be  Irrelevant, 
since  we  will  show  that  colllslonal  time  scales  are  larger  than  the  time  It  takes 
for  an  Ion  to  transit  the  device,  while  magnetic  effects  would  destroy 
convergence  before  a  single  pass.  The  magnetic  deflection  can  be  estimated  by 
calculating  the  orbits  of  an  ion  in  a  magnetic  field  B  -  B0(r/RB^m’  where  rb  *s 
radial  location  where  B  -  BQ,  and  a  radial  potential  $(r),  taking  the  direction 
of  B  to  be  perpendicular  to  14  In  order  to  find  the  maximum  deflection.  This 
was  done  in  Reference  2,  leading  to  a  spatial  spread  In  the  ion  focus  due  to 
magnetic  deflection  from  the  fields  at  r  <  rQ  of 
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where  rQ  is  the  ion  birthpolnt  (rQ  <  Rg)  and  40  the  potential  drop  from  r  »  0 
to  r  -  rQ.  For  the  purpose  of  calculating  collisional  effects,  we  assume  that 
the  source  is  at  an  rQ  which  provides  satisfactory  convergence. 

The  test  Ion  radial  energy  profile  Is  then  approximately 


(3) 
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The  Ion  azimuthal  energy  profile  follows  from  conservation  of  angular  momentum, 
rv,(r)  *  constant,  leading  to 

«r  1  m  J  ,R,2 

Ex  "  ?  V.LR 

where  v±R  Is  the  Ion's  azimuthal  velocity  at  birth. 

We  assume  that  the  Ion  source  Is  active,  for  a  time  long  compared  with  the 
time  it  takes  for  an  ion  to  flow  from  r  «  R  to  r  =  0  and  back  again,  and  that 
the  source  produces  a  low  energy  ion  cloud  at  r  ■  R.  The  ion  distribution 
everywhere,  and  the  radial  Ion  density  profile,  are  then  determined  by  the 
potential  0(r)  and  the  Ion  distribution  function  f^  at  the  birth  point,  which 
can  be  written  in  terms  of  the  constants  of  the  motion  (E  and  rv±),  since  the 
Vlasov  equation  for  the  steady  state  Is  f^  (E,  rv±)  ■  constant.  There  are  three 
quite  different  regions  of  the  device,  characterized  by  the  relative  values  of 
the  Ion's  birth  energy,  and  azimuthal  and  radial  energies.  Within  a  very  small 
radial  distance  Inward  from  the  ion  birth  point,  the  potential  change  is 
negligible,  the  ion  energy  is  of  the  order  of  it's  birth  energy,  and  the  density 
is  given  by  n  -  nedge.  The  potential  Increases  Inward  radially,  and  since  the 
maximum  well  depth  is  in  the  keV  range,  compared  with  birth  energies  in  the 
several  eV  range,  it  follows  that  2e0  »  EQ  a  short  radial  distance  from  the 
edge,  while  n^vjp  (R/r)2  <  2etf.  These  conditions  hold  for  a  large  radial  region 
of  the  device  which  we  refer  to  as  the  bulk.  Evd^ually,  however,  at  a  smaller 
radius,  r  «  R,  m^v^  (R/r£)2  -  2e0max  and  the  ions  can  come  no  closer  to  r  » 
0.  This  defines  the  dense  core  region.  This  inward  motion  clearly  evolves  a 
Maxwellian  source  at  r  -  R  into  a  radially  counterstreaming  bi -Maxwellian  f^  at 
r  <  R  proportional  to 
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where  the  radial  drift  responds  to  the  radial  electric  field,  vd  ■  (2etf/m. J1/2, 
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and  the  transverse  "temperature"  T±  •  m^vxo  (R  /r  )  increases  inwardly  to  provide 
conservation  of  angular  momentum,  with  m^vXQ  an  average  perpendicular  energy  of 
ions  born  at  r  -  R.  If  the  ion  source  is  distributed  rather  than  sharply 
localized,  the  same  arguments  presented  here  can  be  used  with  reference  to 
azimuthal  "slices"  of  the  ion  source. 

Integrating  f^  over  velocity  yields  the  density  in  the  various  radial 
regions  of  the  device.  As  a  model  for  the  ion  distribution  at  the  birth  point, 
we  assume  that  at  birth  the  ion  distribution  function  is  uniform  in  energy  up 
to  some  small  energy  EQ,  and  uniform  in  angular  momentum  up  to  some  small 
azimuthal  velocity  v±(),  and  take  the  potential  at  the  birth  point  to  be  0  -  0. 
This  distribution  is  described  In  terms  of  the  constants  of  the  motion  (E  » 
(l/2)m..(v‘  +  v^)  -  etf,  rv^)  by  the  function 

f i  ■  3<VE0>3/2  "edge  ■  0  <  E  <  Eo  -  0  <  r\  <  R  % 

-  0  else 

A  Maxwellian  gives  similar  results  for  n(r),  but  tijlflat  distribution  is  easier 
to  work  with  analytically.  The  Vlasov  equation  for  steady  state  is  f(E,rvx)  - 
constant,  so  the  density  and  distribution  everywhere  Is  determined  by  the 
potential  and  the  value  of  f  at  the  ion  birth  point.  There  are  three  regions. 
For  a  very  small  distance  from  the  birth  point,  the  potential  is  negligible,  and 
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the  density  is  given  by  n  -  ne(jge.  In  the  bulk  of  the  plasma,  interior  from  the 
edge,  where  EQ  «  2e$,  and  m-jV‘0(R/r)  <  2e0,  the  density  is  given  by  Integrating 
f  over  vxdvxdvr  with  the  result  (rc  is  defined  in  Eq.  (5)  below) 

3  Eo  l'Z  R2 

nB  *  5  <1^  "edge  *  rc  <  r  <  R  (4) 

Assuming  that  the  potential  reaches  its  full  value  at  a  moderately  large 
distance  from  r  *  0,  we  see  that  a  mean  electron  density  throughout  the  bulk  of 
the  device  will  go  as  1/r  .  Moving  inward,  the  (R/r)  factor  eventually  becomes 
substantial,  so  that  the  density  becomes  much  larger  than  nfi(j  .  Eventually, 
the  radius  is  so  small  that  m.v^R/r)  >  2e^max-  Inside  that  radius  the  density 
is  changing  fairly  slowly,  on  the  scale  of  the  changing  potential.  This  defines 
the  radius  of  the  dense  core  rc  as  the  radius  at  which  2e0  -  mivx0(R/r)2’  because 
outside  that  radius  n  -  1/r2,  while  Inside  that  radius  n  is  nearly  constant, 
found  by  integrating  f^  over  velocity, 

Eo  1/2  R  2 

nc  “  Wmax*  ^  "edge  *  r  <  rc 

rc  ■  <V?o/2eW1/2  R  •  <5> 

+ 

The  size  of  the  core  and  the  central  density  are  seen  to  depend  on  the  angular 
momentum  at  the  outer  turning  points  of  the  ion  motion.  This  leads  to  an  overall 
view  of  the  SCIF  which  determines  the  density  and  energy  profiles  needed  to 
calculate  the  collisional  relaxation  of  f^  toward  Maxwellian. 
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B. 


COLLISION  FREQUENCY  PROFILE 


Using  the  ion  distribution  function  discussed  in  the  last  section,  where 
most  of  the  ions  have  turning  points  vr  -  0  at  the  same  radius  r  -  R  at  the 
periphery  of  the  central  potential  ${r),  and  taking  the  radial  and  azimuthal 
energy  moments  of  the  distribution  in  the  three  regions  of  the  device,  we  find 
that  the  core  and  the  edge  are  isotropic,  while  the  bulk  of  the  plasma  is  highly 
anisotropic,  with  Er  »  E  . 

In  the  core,  defined  by  Eq.  (5),  the  ion  density  is  relatively  flat;  using 
the  above  estimates  for  the  ion  energy  and  density,  we  write  the  ion-ion 
(hydrogen)  collision  time  in  the  core,  r^c,  as  (energies  are  in  eV,  density  in 
cm”3) 


2  x  107  (ed  )3/2 
,  v  ymax' 

Tiic  “  n  lnA_ 
c  c 


(6) 


where  the  Coulomb  logarithm  in  the  core  is 


lnAc  =7+2.3  log 
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3/2 


(nc/1014  cm'3) 1/2 


Equation  (6)  gives  r^c  »  1  s  and  30  x  10"6  s,  respectively,  for  the  EXP  and 
REACTOR  parameters  in  Table  I.  In  addition,  lnAc  -  23  and  20  for  the  EXP  and 
REACTOR  parameters. 
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In  the  bulk,  the  Ion  density  decreases  radially  outward  basically  as  1/r  , 
with  an  additional  weak  radial  dependence  from  d,  which  tends  to  make  the  density 
distribution  flatter  than  1/r*  beyond  about  the  R/2  point  of  the  device.  The 
Ion-Ion  collision  time  In  the  bulk,  r^g,  is 


and  using  Eqs.  (4)  and  (5),  which  give 


where 


For  r/R  -  1/2,  InAg  «  27  and  25  for  the  EXP  and  REACTOR  parameters,  respectively. 

The  edge  region  is  defined  as  that  region  in  which  the  ion  radial  velocity 
slows  down  to  a  value  of  the  order  of  it's  birth  speed,  that  is,  Ef  -  e0max  U  - 
(re/R)p]  -  EQ.  The  radial  extent  of  the  edge  region,  Are  -  R  -  rg,  is  then  given 
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by  Arg  «  (l/p)(E</e0nax)R'  The  Ion-Ion  collision  time  in  the  edge,  Tiie,  is 
given  by 
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and  using  Eq.  (5),  this  becomes 
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and  lnAe  -  14  and  7  for  EXP  and  REACTOR  parameters,  respectively. 

The  collision  time  r  Is  not  sufficient  for  comparing  the  effectiveness 
of  collisions  in  the  three  spatial  regions,  because  the  time  an  ion  spends  in 
each  region  is  drastically  different.  A  better  figure  of  merit  is  the  ratio  of 
the  ion  transit  time  ttf  through  the  region  Ar  t^he  collision  time  in  the 
region.  This  essentially  gives  the  number  of  collisions,  Ncoll,  that  a  test  ion 
undergoes  in  a  particular  region  of  the  device  per  transit. 

The  number  of  collisions  in  the  core  per  pass  is  then,  since  rc  - 

V6W1/2r- 
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5  x  10'8  EXP 


=  3  x  10'4,  REACTOR  . 


Thus  one  collisional  upscatter  time  corresponds  to  a  large  number  of  ion 
transits. 

The  ratio  of  the  number  of  collisions  in  the  edge  region  to  the  number 
of  collisions  in  the  core  per  transit  is 


^Ncoll^edge  Are 


M, 


max 


m  TSi 


ic 


Tiie 


which  is  (take  p  *  3) 


^coll^dqe  1  ,e*max,  ,lnAe, 

"  p  ‘tt’ 


4  x  1<T,  EXP 


=  2  x  103,  REACTOR  , 

which  corresponds  to  a  small  fraction  of  an  edge  collision  per  transit  for 
experimental  parameters,  but  to  approximately  one  collision  per  transit  for 
reactor  parameters.  & 

As  representative  values  for  the  bulk  region,  we  take  r  *  R/2.  Then  the 
ratio  of  the  number  of  collisions  in  the  bulk  region  to  the  number  of  collisions 
in  the  core  per  pass  is  approximately 
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|Nco11 |b  _  R/2  T1ic 
(Ncoll>c  rc  Ti  i  B 


» 


which  becomes 


^coll^B  2R  /Jo_w1nAB» 
<Wc  ~rc  (®W(lSV 


5  x  10'2,  EXP 


*  2  x  10'2,  REACTOR  . 


From  the  above,  we  see  that  the  number  of  collisions  per  pass  is  by  far 
the  highest  at  the  edge  region,  where  the  low  energy  ion  distribution  is 
relatively  isotropic.  The  number  of  collisions  is  lowest  in  the  bulk  of  the 
device,  where  the  radial  ion  energy  is  high,  and  where  the  energy  distribution 
is  highly  anisotropic,  with  Er  »  E±. 


III.  HEURISTIC  ESTIMATES  OF  ION  "LOSS"  RATES 


In  this  section  we  give  heuristic  estimates  of  ion  "loss"  rates,  in  the 
sense  of  loss  of  Ion  convergence,  In  a  SCIF  device,  due  to  the  following  ion- 
ion  collisional  effects:  (1)  upscatter  In  ion  radial  energy  in  the  core  region, 
and  (2)  perpendicular  deflection  In  the  bulk  region. 

The  change  in  ion  velocity  by  ion-ion  scattering  in  a  sir.gle  transit 
through  the  bulk  and  core  is  expected  to  be  small,  since  the  ratio  of  transit 
time  to  collision  time  in  the  bulk  and  core  of  the  device  is  small  for  both 
experimental  and  reactor  parameters.  For  example,  from  Section  II,  the  ratio 
of  the  transit  time  to  the  ion-ion  collision  time  in  the  core  is  of  the  order 
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of  5  x  10  for  experimental  parameters,  and  of  the  order  of  -  3  x  10  for 
reactor  parameters.  The  ratio  of  the  transit  time  to  the  ion-ion  collision  time 
in  the  bulk,  for  r  ■  R/2,  is  of  the  order  of  -  3  x  10’9  for  experimental 
parameters,  and  of  the  order  of  -  6  x  10'®  for  reactor  parameters.  Therefore 
we  treat  upscatter  as  a  diffusion  process  (in  velocity  space)  in  the  core  and 
bulk,  rather  than  as  a  single  scattering  process.  In  analogy  with  the  spatial 
diffusion  problem,  we  use  a  continuity  equation  in  velocity  space  as  a  basis  for 
the  loss  process 


dn 

3t 


a2n 

av2 


where  the  diffusion  coefficient  0  -  (6v)  /t,  with  6v  being  the  change  in  velocity 
produced  in  a  time  r  due  to  scattering  in  the  dense  core,  and  n  is  the  time 
dependent  density  of  ions  in  velocity  space. 
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For  the  estimates  In  this  section,  we  assume  that  the  background  Ions  are 
Isotropic  and  are  described  by  a  Maxwellian  distribution;  we  expect  this  to  give 
a  reasonable  result  at  least  to  order  of  magnitude.  The  diffusion  tensor  in  the 
Fokker-Planck  equation,  ■  <Av^AVj>/2r,  Is  diagonal  In  this  case,  and  can  be 
expressed  In  terms  of  components  parallel,  D|,  and  perpendicular,  D±,  to  the 
direction  of  the  test  particle  Initial  velocity  (before  the  collision).  These 
components  are5 


.  eW  ¥^'> 

II  4ire2  m2  v  2(v/v')2 

D  -  eVlnAl_ 

1  8*£0  m2  v 


*i(v/v') 

*  (V/V')  -  - * 

2(v/v'r 


(7) 


where  the  unprimed,  primed  quantities  refer  to  the  test,  field  ions, 
respectively;  the  temperature  T'  for  the  field  ions  Is  (3/2)  T'  *  (1/2)  n^v'2, 
InA'  is  the  Coulomb  logarithm  for  the  field  ions,  and  §j(x)  is  defined  in  terms 
of  the  error  function 


*(x)  -  ^  J  exp  (-f2)  df 
Jx  0 


by 


*l(x) 


♦(x)  -  —  exp(-x2)  -  *(x)  -  x  || 

Jx 


For  x 


1,  *j/2x2  -  0.214,  and  *  -  0.843. 6  For  x  >  2,  fj  «  #  «  l.5 
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A. 


ION  LOSS  DUE  TO  RADIAL  UPSCATTER  IN  THE  CORE 


Upscatter  In  radial  velocity  In  the  core  by  collisions  can  affect  ion 
convergence  In  the  device.  Consider  a  collision  which  causes  one  ion  to  gain 
energy  and  another  Ion  to  lose  energy.  The  ion  which  gains  energy  in  the  core 
can  make  a  larger  radial  excursion  to  a  radius  beyond  it's  birth  radius,  with 
the  possibility  that  the  magnetic  fields  at  the  new  turning  point  will  be 
sufficient  to  deflect  the  ion  in  v^,  leading  to  a  degradation  of  the  focus. 

We  describe  radial  upscatter  in  the  core  by  a  diffusion  equation 


dn 

at 


Here  D|e^  -  D|c  (tc/tD),  where  D|c  is  the  parallel  diffusion  coefficient  in  the 
core,  given  by  Eq.  (7)  with  1/2  n^v'2  «  etfmax  and  n'  «  nc,  and  (tc/tp)  is  the 
fractional  time  that  a  test  ion  spends  in  the  core  per  transit  through  the 
device.  Assuming  the  core  has  radius  r  ,  with  uniform  density  n  ,  the  transit 

V#  L 

time  through  the  core  Is  given  approximately  by  tc  ■  2rc/vc>  where  vc  is  the 
velocity  of  the  test  Ion  in  the  core  region.  The  time  for  a  test  ion  to  transit 
through  the  entire  device  is  given  approximately  by  tp  -  2rQ/vr,  where  rQ  is  the 
radial  birth  location  of  the  test  Ion,  and  vr  Is  the  average  velocity  of  the  test 
ion  in  the  device,  given  by  vf  -  (1/R)  J  (Ze^max/|lp1^2  U  •  (r/R)p) i/<2  dr  -  x/4 

(2e^max/mi)1/2  ~  */4  V  for  P  *  2-  Thus 
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where 


T|c  -  2«JmJv'3/(e4nclnAc) 

Since  the  test  particle  Is  actually  one  of  the  core  particles,  we  assume 
that  v  -  v'(l  +  f),  where  f  -  Av/v'  and  Av  is  the  upscatter  in  the  test  particle 
velocity  due  to  collisions  in  the  core.  Then 


[#j(l  +  f)] 


1 


(1  +  f)’ 


r  v  n»(1012cm  3) 
(^)<r)  (10  kV) 


We  estimate  the  ion  "loss"  rate  rlos$  due  to  upscatter  in  radial  velocity 

using 


(Avr)2 

Ti°s$  ■  &|efr 

where  Avr  *  fv;  is  the  increase  in  the  test  ion  radial  velocity  due  to  collisions 
In  the  core,  which  we  relate  to  the  loss  of  focus  below.  The  loss  time  due  to 
upscatter  would  then  be  given  by 
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s 


.  (8) 
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Since  the  perpendicular  deflection  of  the  Ion  by  the  magnetic  field,  AV, 
at  the  ion  birth  location  Is  related  to  the  core  radius  by2 


AV 

r 


♦ 


the  Increase  In  core  radius,  Arc,  due  to  larger  deflection  at  r^  is  then 

Arc  AV(r^)  -  AV(rQ) 
rc  * 


The  perpendicular  deflection  of  the  Ion  by  the  magnetic  field  in  the  SCIF  device 
was  estimated  In  Ref.  2  by  considering  Ion  orbits  in  orthonormal  magnetic  and 
electric  fields: 


where  $Q  is  the  potential  drop  from  r  ■  0  to  r  ■  rQ.  At  r^, 

AV(r')  «|f  ,  (1  +  f)2(2m+2/p)r2m+2 

~T~  “  <7TT?  r®  ’ 

* 

so  that  the  change  In  AV  due  to  a  larger  radial  excursion  to  is,  for  f  <  1, 
AV(ro>‘  W0)  ^  2 (2m  +  2)f 

rc  “  P 
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Thus  Arc/rc  >  1  requires  f  comparable  to  unity.  We  estimate  the  ion 
"loss"  time  from  Eq.  (8)  evaluated  with  f  -  1.  For  rQ  «  R,  and  for  the  EXP 
parameters  given  in  Table  I, 


TloSS  -  450  s 


For  REACTOR  parameters, 


Tloss 


45  ms 


* 


as  compared  with  discharge  time  scales  of  -  10  ms. 


B.  ION  LOSS  DUE  TO  PERPENDICULAR  DEFLECTION  IN  THE  MANTLE 


Perpendicular  deflection  due  to  scattering  in  the  bulk  of  the  device  could 
lead  to  isotropization  of  the  anisotropic  ion  distribution.  From  conservation 
of  angular  momentum  (without  scattering),  rv^(r)  *  constant.  An  increase  in  the 
azimuthal  ion  velocity,  Av±,  can  be  related  to  an  increase  in  the  core 
convergence  radius,  Arc,  by2 


Avx(r) 

vjtrT 


Thus  isotropization  in  the  bulk  region,  resulting  in  a  transfer  of  vr  -  v  in 
the  bulk,  could  degrade  the  Ion  focus.  In  this  section,  we  estimate  the  rate 
at  which  Ions  are  "lost"  to  the  dense  core  in  the  sense  that  they  converge  to 
a  significantly  larger  core  radius. 
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We  use  a  continuity  equation  In  velocity  space  as  the  basis  for  the  loss 


process 


dn 

3t 


i.eff 


A 

5? 


> 


where  is  the  effective  perpendicular  diffusion  coefficient  in  the  mantle 


Dieff  “  Di  {W 


w#ith  D±  the  perpendicular  diffusion  coefficient  in  the  bulk,  given  by  Eq.  (7) 
with  T#  -  etfmax  and  n'  -  ng,  and  (tg/tp)  is  the  fractional  time  that  a  test  ion 
spends  in  the  bulk  per  transit  through  the  device. 

As  representative  values  for  calculating  the  loss  rate  due  to  deflection 
in  the  bulk  of  the  device,  we  use  r  «  R/2,  and  we  assume  the  background  ions  in 
the  bulk  region  have  a  uniform  density  ng  (r  ■  R/2)  given  by  (see  Eq.  (4)) 


nB  (r  -  R/2) 


f  E. 


3n 


edge 


ed 


max 


1/2 


» 


and  a  radial  energy  given  by  Er  -  «0max*  The  transit  time  for  a  test  ion  in  the 
bulk  region  Is  then  of  order  of  1/2  its  transit  ^me  through  the  entire  device. 
The  effective  perpendicular  diffusion  coefficient  is  then 


D 


1  JL 

ieff  ~  5  T|B 


* 


»1(v/v,)~ 

2(v/v')2 


* 
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where  T|b  -  tjc  (n^ngHlnA^lnAg)  -  1/3  (e$max/E0)  r|c,  and  v'  »  v'  = 

(2e0max/mi)  ^  . 

We  estimate  the  ion  "loss"  time,  in  the  sense  of  ion  focus  degradation, 
as  the  time  for  (Av  /v  ),  and  thus  (Ar  /r  ),  to  become  of  order  q,  where  q  is 
a  number  comparable  to  or  greater  than  1.  By  conservation  of  angular  momentum, 
v±  -  v±0  (R/r)  in  the  bulk,  so  that  Lv±  -  qv±  -  2qv1Q  at  r  -  R/2.  The  ion  loss 
time  can  then  be  estimated  from 


loss 


(Av^ 

Dieff 


-  11 


q2r 


(Jr) 


1 


*  (St)  - 


*j(v/v')" 


2(v/v')4 


18  q2  r 


for  v  ■  v' 


(9) 


with 


rjc  -  0.3  sec 


t^ax<10  keV))3/Z 
nc(1012cm'3} 


For  the  experimental  parameters  in  Table  I,  the  loss  time  is  of  order  5 
q  sec.  For  reactor  parameters,  the  ion  loss  time  due  to  this  collisional 
process  Is  of  the  order  of  0 

Tloss  -  02  q2  “ 

as  compared  with  discharge  time  scales  of  -  10  ms. 
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2 

Note  that  the  loss  time  Is  proportional  to  the  (Avx)  required  for  loss, 

o 

which  Is  proportional  to  the  (Arc)  required  for  loss.  The  value  of  Avx  required 
for  loss  will  be  calculated  later  In  this  report. 

C.  THERMAL IZATION  BY  COLLISIONS  AT  THE  EDGE 


A  test  Ion  spends  approximately  the  same  amount  of  time  transiting  the 
core  region  as  it  does  transiting  the  thinner  edge  region.  That  is,  the  transit 
time  In  the  core,  tc,  where 


tc  -  2rc/vc  -  2R 


K  1 

1/2 

■i  1 

ed> 

L  ymaxJ 

1/2 
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Is  comparable  to  the  transit  time  In  the  edge  region,  te(j  ,  where 


edge 


2Are  2R 

- > 
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o 
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vo  P 
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1/2 
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From  Section  II,  the  Ion-ion  collision  frequency  is  larger  at  the  edge  than  in 
the  core,  viz., 


Vie  estimate  whether  the  ion  distribution  might  thermal ize  in  the  edge 
region.  The  thermal Izat Ion  time  scale  Is  of  the  order  of  the  Ion-ion  self 
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collision  time  r^e.  The  amount  of  time  that  a  test  Ion  spends  In  the  edge 
region  Is  tfi^e  times  the  number  of  passes  through  the  device. 

-4 

For  experimental  parameters,  f^e  -  1  sec  •  (Eo/e0max)  -  5  x  10  sec. 

The  Ion  transit  time  through  the  edge  region  would  be  of  the  order  of  10  ns  for 

R  -  100  cm.  Assuming  10*  transits  through  the  device,  a  test  ion  would  spend 

about  10'*  sec  in  the  edge  region,  somewhat  less  than  the  ion-ion  collision  time. 

-5 

For  reactor  parameters,  however,  r-ie  -.(3  x  10  sec)  •  (^0/e^max)  -1*5 
ns.  The  ion  transit  time  through  the  edge  region  would  also  be  of  the  order  of 
1.5  ns  for  R  -  200  cm,  p  -  2.  Thus  in  several  passes  through  the  device,  the 
ion  distribution  at  the  edge  could  be  thermal ized. 


IV.  CALCULATION  OF  THE  ION  DISTRIBUTION  PRODUCED  BY  COLLISIONS 


In  Section  II  we  calculated  the  relative  rates  of  collisions  in  the 
various  spatial  regions  of  a  SCIF,  and  found  that  collisions  in  the  edge  regions 
dominate.  More  importantly,  we  found  that  the  collision  rate  in  the  interior 
of  the  device  is  very  small  compared  with  the  transit  frequency,  so  that  for  many 
transits  the  ions  return  to  their  original  turning  points  at  r  -  R.  That  is, 
the  definition  of  "edge  region"  is  preserved  for  thousands  of  interior  plasma 
collision  times.  In  Section  III  we  estimated  collisional  ion  "loss"  time  scales, 
and  found  that  these  time  scales  are  much  longer  than  discharge  time  scales  for 
experimental  parameters,  but  not  for  reactor  parameters.  We  also  found,  however, 
that  collisions  in  the  edge  region  might  relax  the  anisotropy  in  that  region. 
In  this  section  we  give  a  Fokker-Planck-like  calculation  of  temperature 
relaxation  of  the  edge  distribution  due  to  ion-ion  collisions  throughout  the 
device.  Since  the  collision  frequencies  In  the  bulk  and  core  are  much  lower  than 
the  transit  frequency,  the  edge  distribution  determines  the  evolution  of  the  ion 
density  and  energy  structure  in  all  three  regions.  The  edge  distribution  evolves 
(a)  due  to  collisions  at  the  edge,  (b)  due  to  the  influx  of  energy-scattered  ions 
(due  to  collisions  in  the  core)  and  (c)  due  to  energy-deflected  ions  (due  to 
collisions  In  the  bulk)  coming  into  the  edge  region  in  every  transit. 

The  dominant  effect  of  edge  collisions  ^  to  relax  any  temperature 
anisotropy  which  develops  in  this  region.  The  source  of  ion  anisotropy  in  the 
edge  is  the  influx  of  excess  v  into  the  edge  region  due  to  the  conversion  of 
radial  to  azimuthal  energy  via  ion-ion  collisions  in  the  highly  anisotropic  bulk 
plasma.  The  goal  of  the  present  calculation  is  to  quantify  the  extent  to  which 
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Isotropizlng  collisions  In  the  low  Ion  energy  edge  can  remove  or  reduce  this 
Increase  of  v  .  The  anisotropy  at  the  edge  Is  the  primary  effect  which  spreads 
out  the  distance  of  closest  approach  to  r  -  0  (this  is  obvious  from  energy 
conservation),  which  In  turn  could  lead  to  an  isotropic  interior. 


A.  ION  DISTRIBUTIONS  IN  THE  THREE  REGIONS 


In  Sections  IV. A  and  IV. B,  we  assume  that  the  effect  of  the  competition 
between  edge  collisions  and  bulk  collisions  will  be  that  the  edge  region  will 
ultimately  tend  toward  a  bi -Maxwell Ian,  characterized  by  Tf  and  T  ,  and  calculate 
the  steady  state  difference  of  Tr  and  T±  in  the  edge  region. 

Thus  for  the  edge  region,  we  take 


Kee 


Vr 

2T„ 

ro 


_  Ji 

V± 


10 


(10) 


with  Ke  being  a  normalization  constant  such  that 
"edge  *  f  ^  X  fi,e 

Since  the  Ion  collision  frequency  in  the  interior  is  low,  the  ions  maintain  a 
dominant  set  of  turning  points  at  r  -  R,  so  that  the  interior  distribution 
continues  to  be  the  evolution  of  an  edge  Maxwellian  under  the  influence  of  a 
strong  central  potential  0(r).  Therefore,  we  take  the  ion  distribution  in  the 
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bulk  region  to  be  the  sum  of  two  radially  counter-drifting  bl -Maxwell Ians.  The 
radial  drift,  Vg  «  /2e^(r)/m^ ,  Is  due  to  the  electrostatic  field.  The  azimuthal 
temperature  T  Increases  with  decreasing  radius,  reflecting  conservation  of 
angular  momentum,  T±(r)  -  Txo  (R/r)  .  Thus  in  the  bulk 


(11) 


where  Kg  is  a  normalization  constant  such  that 

nB--fd3Xf1.B  • 


We  note  that  the  bulk  distribution  reduces  to  the  edge  distribution  at  r  *  R, 
since  0{R)  -  0. 

We  assume  that  the  core  is  defined  as  that  region  in  which  the  azimuthal 
and  radial  ion  energies  are  the  same,  that  is,  Tf  *  T±  «  e^max'  Thus  we  mode1 
the  core  distribution  as  an  isotropic  Maxwellian  at  the  temperature  T£  «  e0max> 
since  0  is  virtually  constant  over  the  small  core  region, 


m^v2 

'  2TT 

f1,c  ■  Kce  #  <12> 

For  simplicity,  we  use  the  density  and  energy  profiles  in  the  three 
regions  of  the  device  derived  in  Section  II,  which  can  be  shown  to  be  similar 
to  those  produced  by  a  bi -Maxwellian. 
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B. 


CALCULATION  OF  THE  STEADY-STATE  ION  DISTRIBUTION 


We  model  the  evolution  of  T  In  the  edge  region  by: 


dTio  .  dTio 
dt  dt 


+  dTio 
+  TT 


ISO 


(13) 


Here,  dTXQ/dt | g  Is  the  source  term  for  the  additional  v  coming  in  to  the  edge 
region  on  each  ion  transit,  due  to  the  conversion  of  radial  to  azimuthal  energy 
by  collisions  in  the  bulk,  and  dTxo/dt | j^q  is  the  isotropization  rate  due  to 
collisions  in  the  edge  region.  There  is  a  similar  equation  for  T  .  However, 

2  2  7 

since  (mv£/2)  -  Tp/2,  m v^/2  -  T  ,  we  can  use  conservation  of  energy  to  write 


dT 

~ar 


d  ^ro 

•  3tT 


(14) 


and  to  obtain  an  equation  for  the  evolution  of  the  edge  temperature  difference, 


AT 


dT 


dT 


+  3 


10 


"3T 


ISO 


(15) 


There  is  a  similar  driver  In  which  Tro  is  changed  by  energy  scattering 
in  the  core,  with  Txo  in  the  edge  changed  becau#  energy  conservation  slightly 
alters  the  turning  point.  We  emphasize  here  the  scattering  because  T± 
acquired  during  a  transit  is  carried  out  to  the  turning  point  decreased  only  by 
r/R  (of  order  1/2),  while  vf  always  goes  to  0  at  the  turning  point  because  of 
the  structure  of  0(r). 
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r/R  (of  order  1/2),  while  vr  always  goes  to  0  at  the  turning  point  because  of 
the  structure  of  tf(r). 

In  the  following  calculation,  we  assume  that  the  edge  temperature 
difference  is  small,  AT  -  T  -  Tro  «  T,  because  the  conversion  of  radial  to 
azimuthal  velocity  of  an  ion  in  each  transit  through  the  bulk  is  small  compared 
with  the  rate  of  isotropizing  edge  collisions,  from  the  estimates  in  Section  II. 

The  isotropization  time  for  a  bi -Maxwellian  distribution,  in  the  limit 
of  small  temperature  differences,  is  given  by:7 
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(16) 


where  rg  «  (15/4)  (m{/2  T^//*  e4  nedge  lnAe),  and  Tfidge  -  Tro  to  this  order 
of  approximation. 

To  calculate  (dT10/dt)  due  to  the  bulk  region,  where  the  ion  energy 
distribution  is  highly  anisotropic,  we  write  the  Fokker-Planck  collision  term 
using  the  Landau  form  for  the  collision  operator.  We  take  the  energy  moment  of 
the  kinetic  equation 


at  “  5v7 
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where  j  is  the  three-dimensional  particle  fl#  density  in  velocity  space, 
representing  particle  collisions: 
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where 


■ft  r 


<AV  ‘  <AvftAV 


dafi  “  h  <AvftAV 


as  given  for  example  in  Reference  8. 


Using  conservation  of  energy,  we  obtain  for  the  increase  of  TiQ  due  to 
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We  compute  jr  by  writing  the  flux  in  this  case  in  Landau  form9 
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with 
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Here  the  primed,  imprinted  quantities  refer  to  the  field,  test  particles, 
respectively,  and  u  -  |v  -  v'|.  We  assume  that  the  field  and  test  particles 
counter-drift  radially: 
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which  can  be  written  as 
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and,  using 


2k  f  dv  J 

-0  0 


9 


Eq.  (19)  becomes 
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dvrdvr 
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where  A  Is  given  by 
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m1  ,2  InA  /4jre2.2 
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To  do  this  Integral,  we  make  the  approximation  that  u  -  I u I  -  |v  -  v'l  -  2  v„. 

i  —  i  *  r  r 1  D  ’ 

so  that  u  can  be  taken  out  of  the  Integral.  We  rewrite  the  Integration  variables 
as 


*; 


v;  +  vd 


9 


v  -  v ' 
r  r 


*r  "  *r  +  2*0 


to  obtain 
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~3T 
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ro 


miui 

2Vn)  -  4^  U. 


"i*5 


m  u '  2 

Vr 

2T~ 

ro 


Vj. 

*7 


dvr  dvj.  dv2  dvj2,  (23) 


which  can  be  rewritten  as 
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In  the  integrations  over  dvr,  dv^.,  only  those  terms  which  are  proportional  to 
even  powers  of  vr,  v^.,  are  non-zero;  thus  the  integral  over  dvr,  dvj.  is 


*!  “  I  dvr  /  dv^ 

-to  •to 


miu?  -2  miu?  -2  2 

t1^  v2  -  (v2  +  2v2) 
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which  is 
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so  that 
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Tro  2m1vD 


m1vi  m1vi 


/2 


2T~  '  TT7  -  2 
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2  2 

To  do  the  v^,  Integrations,  we  write  *  (v±  -  v^)  ,  so  that  the  integration 
over  dv±,  dv^  becomes 


In  - 


*  2  " 


f  dvz  f  dv'2  (v2 

0  1  o  x  1 
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and  we  have 
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Then  with  the  normalization  constant  Kg  given  by 


nij 

kb  ‘  nB  (27r)  J 

i  ro 


we  have  that 


akr  2*T  2T,  3  ,  nR  2jtT  - 

3  <-ir^  <isr>  <2  -  ?>  -  ~2  7ZT3 lnA  e  <4  -  *) 


(2Vn) 


2  m2  (2VD)' 


We  rewrite  Eq.  (24)  as 


dT±o  (2,Vc  T  Tro  ‘  V 
“ar  "  rB 


(25) 


where  the  Isotroplzation  time  xB  Is  given  by 


m2  (2v0)3 


B  2*e4  (InAg)  ng  (4  -  *) 


(26) 


In  order  to  use  this  as  a  source  term  for  (dT/dt)g  in  the  edge  region. 
It  Is  necessary  to  take  Into  account  the  fact  th^an  ion  in  it's  transit  through 
the  device  spends  more  time  In  the  bulk  than  In  the  edge.  Thus,  we  multiply  Eq. 
(25)  by  the  ratio  of  the  transit  time  in  the  bulk  to  the  transit  time  ax  the  edge 
(tg/te).  This  is  equivalent  to  writing  Eq.  (25)  In  units  of  the  edge  transit 
time.  We  also  need  to  take  the  following  effect  into  account,  which  is  a 
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consequence  of  conservation  of  angular  momentum.  An  ion  which  gains  excess 
2 

(l/2)m^v2  in  the  bulk  at  radius  r,  will  bring  an  excess  azimuthal  energy  into 

2  2 

the  edge  of  magnitude  *  (l/2)m^v*(r/R)  .  We  thus  have,  for  use  in  Eq.  (13), 


dT.o 

~ar 


.  2"lv0  ,r,2  H 


(27) 


Equation  (15)  then  becomes 
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3t 


AT 


6Vd 

tB 


(28) 


We  use  the  estimates  for  the  transit  time  in  the  bulk  and  the  edge  regions  from 
Section  II.  The  ratio  of  the  transit  time  in  the  bulk  to  the  transit  time  in 
the  edge  Is  of  the  order  of  tB/te  -  W(EQ/edmax)(R/p)]  t(Ej/2)/(edmax)1/2]  - 
p/e^max^ro’  w^e  the  ratio  of  the  isotropization  time  in  the  bulk,  rg,  to  the 
isotropization  time  in  the  edge,  re,  is 


f-!/2  (^x>3/2  1 

'4  "edge  1nAe 
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2ne  ng(4  -  *)lnAg 
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which  becomes 
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edge 


(29) 


38 


Thus,  as  long  as  n^gg/ng  >  P»  with  P  defined  in  Eq.  (2)  typically  of  order  2-3, 

we  expect  that  the  temperature  difference  in  the  edge  region  would  remain  small, 

AT  <  TrQ.  This  condition  depends  on  the  form  of  the  density  profile  in  the 

device,  which  Is  determined  from  the  edge  Ion  distribution  and  conservation  of 

energy  and  angular  momentum.  Using  the  flattop  ion  distribution  discussed  in 

Section  II,  we  have  that  nedge/nB  -  (e$max/E0)1/2  »  1. 

The  equilibrium  solution  of  Eq.  (28)  yields  AT  -  (T±  -  Tr)edge  -  Tro 

<Tro/e*max  «  TrQ.  Thus  the  quasisteady  solution  to  this  nonlocal  collision 

relaxation  process  is  an  ion  distribution  which  is  nearly  Maxwellian  at  the  edge. 

Since  the  ions  have  turning  points  in  the  same  vicinity  as  their  birth  point, 

the  interior  distribution  continues  to  evolve  from  the  edge  distribution, 

conserving  energy  and  angular  momentum,  so  that  the  non-Maxwell ian  distributions 

2 

Eqs.  (21)-(22)  with  T±  -  T10  (R/r)  continue  to  describe  the  bulk  plasma.  The 
interior,  like  the  edge,  is  driven  toward  isotropy  by  local  collisions,  but 
driven  toward  anisotropy  by  influx  of  ions  from  the  Maxwellian  edge.  It  is 
because  the  edge  collision  frequency  per  transit  is  many  orders  of  magnitude 
larger  than  the  bulk  rate  that  the  anisotropic  influx  so  dominates  the  interior. 


V.  APPLICATION  TO  SCIF  ION  LOSS  RATES 


For  the  experimental  parameters  in  Table  I,  the  rate  of  collisions  per 
transit  is  so  small  as  to  be  negligible.  For  reactor  parameters,  however,  the 
collision  rates  must  be  compared  to  the  fusion  energy  multiplication  rate,  since 
about  10  transits  before  loss  are  required  for  economic  operation.  Therefore 
in  this  section  we  consider  in  detail  the  loss  produced  by  collisions  for  reactor 
parameters  in  a  magnetically  confined  SCIF. 

The  previous  results  show  that  edge  collisions  maintain  a  non-Maxwellian 
interior  for  a  time  larger  than  the  small  angle  scattering  time  which  in  most 
cases  produces  a  local  Maxwellian.  Applying  this  result  to  SCIF  loss  rates,  we 
note  that  there  are  two  effects.  First,  large  angle  scattering  will  prevent  ions 
from  returning  to  the  vicinity  of  their  birth  location,  thus  removing  the  healing 
effect  of  edge  collisions  which  assumes  a  fixed  range  of  turning  points.  Second, 
the  slightly  anisotropic  edge  distribution  found  In  Section  III  will  produce  some 
degradation  in  the  core  on  the  time  scale  of  the  small  angle  scattering  time  in 
the  core. 

A-  ION  LOSS  BY  LARGE  ANGLE  SCATTERING 


The  frequency  of  90°  scattering  is  obtained  by  simply  setting  log  A  -  1 
in  the  collision  rates  In  Section  III.  To  e^mate  the  loss  rate  we  must 
estimate  how  large  a  change  in  velocity  Av  is  required  to  cause  an  ion  to  be  lost 
to  the  core.  In  this  section  we  consider  in  detail  such  losses  for  reactor 
parameters  in  a  magnetically  confined  SCIF. 
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1.  Ion  loss  Due  to  Radial  Upscatter  In  the  Core 


Upscatter  In  radial  velocity  in  the  core  by  collisions  can  affect  ion 
convergence  in  the  device.  Consider  a  collision  which  causes  one  ion  to  gain 
energy  and  another  ion  to  lose  energy.  The  ion  which  gains  energy  in  the  core 
can  make  a  larger  radial  excursion  to  a  radius  beyond  it's  birth  radius,  with 
the  possibility  that  the  magnetic  fields  at  the  new  turning  point  will  be 
sufficient  to  deflect  the  ion  in  v^,  leading  to  a  degradation  of  the  focus. 

We  estimate  the  value  of  6v  which  could  lead  to  ion  "loss,"  in  the  sense 
of  a  broadened  core  radius  rc,  due  to  increased  deflection  by  the  higher  B-field 
as  the  upscattered  ion  makes  a  larger  radial  excursion.  To  do  this,  we  estimate 
the  new  radial  location  r^  at  which  the  upscattered  ion's  radial  energy  decreases 
to  its  birth  energy.  Neglecting  scattering,  we  assume  the  ion  radial  energy 
profile  is,  generalizing  Eq.  (3)  to  rQ  *  R, 
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Eo  + 


W, 


max 


(30) 


Due  to  upscatter  in  the  core,  the  ion  receives  a  kick  in  energy  AEr, 


roP 


AEr  -  e*max  <1T>  f<2  +  f> 


(31) 


where  f  -  Avr/vr. 
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or 


r 

K 


-  (1  +  f)2/p 


(32) 


If  this  new  turning  point  r'Q  takes  the  ion  into  a  region  of  high  magnetic 
field,  B  -  B0(r/Rg)m,  the  perpendicular  deflection  of  the  ion  by  this  magnetic 
field,  Av±,  changes  the  core  radius  by  (see  Eq.  (5)) 


r 


c 


R 


(33) 


The  perpendicular  deflection  of  the  ion  by  the  magnetic  field  can  be  estimated 
from  straightforward  consideration  of  ion  orbits  in  orthonormal  magnetic  and 
electric  fields.  Assuming  that  the  ions  were  originally  born  at  a  location  rQ 
at  which  magnetic  deflection  would  give  a  contribution  to  h\i±Q  (and  to  rc)  of 
the  same  order  of  magnitude  as  the  hv±  arising  from  the  thermal  spread  T^q  in 
the  ion  source,  it  is  easy  to  show  that  the  spread  Av  due  to  a  larger  radial 
excursion  to  r^,  at  which  the  magnetic  field  is  larger  (B  -  rm),  is 


Av<r;>  4m  ^r  e  ^c 
Av(rQ)  p  vr  “  rc 


(34) 


Thus  ion  loss,  In  the  sense  that  Arc/rc  -  1,  requires  6wr/yr  comparable 
to  unity,  which  means  a  single  large  angle  scattering.  From  the  results  in 
Section  II,  and  the  reactor  parameters  in  Table  I,  this  gives  (setting  log  A 

-  1) 
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(35) 


'loss  -  r1ic  1nAc  rc-r\H  1nAc  -  85  ms 

which  Is  large  compared  to  fusion  reaction  times  of  order  10  ms,  from  References 
1  and  2. 


2.  Ion  Loss  Due  to  Perpendicular  Deflection  in  the  Bulk 


Perpendicular  deflection  due  to  scattering  in  the  bulk  of  the  device 
could  lead  to  isotropization  of  the  anisotropic  ion  distribution.  From 
conservation  of  angular  momentum  (without  scattering),  rv^(r)  -  constant.  An 
increase  in  the  azimuthal  ion  velocity  Av±  at  its  birth  point  is  related  to  an 
increase  in  the  core  convergence  radius,  Arc,  by 


i<ro>  4rc 


Since  vf  »  v±  through  much  of  the  ion  orbit,  even  a  small  decrease  in  Avr/vf 
could  produce  enough  Av±/vx  -  (A vr/vr)(vr/v1)  to  degrade  the  ion  focus.  However, 
a  small  Avr  would  still  result  in  a  turning  point  at  r  «  rQ,  where  collisions 
would  restore  f^  as  calculated  in  Section  III.  In  order  to  change  rQ 
significantly,  6vr/vr  «  1  is  required,  even  for  collisions  in  the  bulk.  The 
result  must  take  into  account  the  lower  density  ng/nc  in  the  bulk,  and  the 
extended  dwell  time  R/rc  in  the  bulk.  This  gives  an  estimate  of  the  loss  time 
due  to  collisions  in  the  bulk  of 


Tloss,B  *  TloSs,core("</"B)<VR>  *  Tloss,core(e*/Eo>1/2  “  10  5  •  <37> 
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for  reactor  parameters.  This  Is  negligible  compared  with  *fUSi0n- 
B.  Loss  Due  to  the  Edge  Anisotropy 


In  Section  III  we  calculated  the  anisotropy  which  develops  at  the  edge 
due  to  small  angle  scattering  In  the  rest  of  the  device.  In  this  section  we 
estimate  the  degradation  of  the  core  which  results  from  this  anisotropy. 

From  Section  III  we  showed  that  the  increase  in  T,„  at  r  *  R  was  of  order 
Eo^/e*)1/2,  and  that  the  Increase  in  core  size  was  Arc/rc  <*  Av  /v±.  Therefore 
it  follows  that 


Arc 

r  v 
c 


y/4 < 

i?)  < 


10 


-1 


(38) 


with  this  small  Ar£  taking  place  on  a  time  scale  of  a  few  ion  transits.  We 
conclude  that  small  angle  scattering,  when  the  Maxwellization  of  the  plasma  near 
the  ion  turning  points  is  taken  into  account,  is  a  very  small  contribution  to 
core  degradation,  even  though  this  is  the  fastest  collisional  effect. 
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Table  I.  Representative  SCIF  parameters. 


EXP 

REACTOR 

e0max  (depth  of  the  potential  well) 

10  keV 

100  keV 

Eq  (ion  birth  energy) 

5  eV 

5  eV 

n£  (ion  core  density) 

1012  cm’3 

1018  cm'3 

R  (radius  of  device) 

1  m 

2  m 

1  T 


B  (magnetic  field  at  R,  m  -  3) 


0.2  T 


